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Abstract. We discuss INTEGRAL'S capabilities to detect a high redshift population of Gamma-Ray Bursts (GRBs) in compar- 
ison to other high-energy missions. Emphasis is done on the study of the relative capabilities of IBIS on board INTEGRAL with 
respect to SWIFT and HETE 2 for detecting a high redshift population of GRBs. We conclude that, if the GRB rate is propor- 
tional to the star formation rate, INTEGRAL'S capabilities for studying GRBs are complementary to the ones of missions like 
SWIFT and HETE 2, specially devoted to prompt localisations of GRBs. Whereas SWIFT and HETE 2 would detect a higher 
number of GRBs than INTEGRAL, IBIS might be able to detect high redshift (z > 7) GRBs, unreachable for SWIFT and 
HETE 2. We discuss the relevance of performing near-infrared (NIR) observations of the INTEGRAL GRBs and the strategy 
that large class telescopes might follow. 
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1. Introduction 

Gamma-Ray Bursts (GRBs) appear as brief (10^'^ s < duration < 
1 000 s) flashes of cosmic high-energy photons, emitting the bulk of 
IL* their energy above ~ 0.1 MeV. Among the thousands of GRBs de- 
. ^ tected in 7-rays since 1967 (Klebesadel et al. 1973 ( only ~ 50 GRBs 
have been identified at optical wavelengths'. For 36 of them spec- 
5_j troscopic redshifts have been measured, ranging from z — 0.0085 
(Galamaetal.l99Fl to 2 = 4.50 (Andersen etal.'2000'(. 

A canonical model has emerged for powering long GRBs (dura- 
tions > 2 s): collapse of a massive star onto a black hole. This leads 
to an intense flash of 7-ray photons followed by an expanding fireball 
which emits radiation at lower frequencies (Woosley 1993 1. The spec- 
troscopic association of GRB 030329 with the type L SN 2003dh sup- 
ports strongly this framework (Hjorth et al. 12003 1 Stanek et al. 20Q3 (. 
Short GRBs (durations < 2 s) have not shown yet a conclusive optical 
counterpart (see Castro-Tirado et al.'2002'for a further discussion). 

We know today that most long duration GRBs originate at 2: > 
0.5. In fact, it is thought that their intrinsic brightness would allow 
us to detect these events at epochs corresponding to the formation of 
the earliest stellar populations. Thus, they may be used as probes into 
the first stages of star formation and their spectra may reveal the early 
heavy-element enrichment of the interstellar medium (ISM). GRBs 
originated by the exploding primitive population of stars are suggested 
to be detectable up to redshifts z; ~ 30 (Lamb & Reichart ii2000> . 

ESA's INTEGRAL satellite offers unique capabilities for the de- 
tection of GRBs thanks to its high sensitivity and imaging capabil- 
ities at 7-rays, X-rays and optical frequencies. INTEGRAL is the 
first 7-ray spacecraft that combines imaging instruments of high pre- 



cision and a continuous real time telemetry link. Developments in 
the GRB field over the past few years have made increasingly clear 
that INTEGRAL, and specially the IBIS" 7-ray imager, could make 
a significant contribution. IBIS is the coded mask telescope on board 
INTEGRAL dedicated to imaging over a wide (15 keV —10 MeV) en- 
ergy range (Ubertini et aI. I2003> . To date, neither of the three classical 
INTEGRAL GRBs (considering GRB 031203 as an X-ray flash) with 
optical counterparts identified have had their redshift measured (Gotz 
et al. l2003l Castro-Tirado et al. l2003l Masetti et al. l2004l . 

In Sect.|2|we perform a rough comparison of INTEGRAL'S sen- 
sitivity with respect to other 7-ray and X-ray missions-instruments. 
Sect. |3| studies in more detail the capabilities of the most sensitive 
high-energy missions-instruments (INTEGRAL/IBIS, HETE 2AVXM 
and SWIFT/BAT) to detect high redshift GRBs. Further, in Sect.|4|we 
focus on the INTEGRAL/IBIS and SWIFT/BAT relative number of 
high redshift GRB detections. As an ado, Sect.|5|makes emphasis on 
the relevance of the near-infrared (NIR) observations to detect high 
redshift afterglows. Last, Sect.|S|presents the conclusions of our study. 



2. Sensitivity comparisons: INTEGRAL 
expectations 

In TableQ we provide the main performances of a number of different 
space missions/instruments with capabilities for GRB research. The 
missions/instruments are divided in two groups depending on their en- 
ergy range. In the first group the INTEGRAL/IBIS sensitivity is nor- 
malized with respect to the Burst Alert Telescope (BAT) 7-ray instru- 
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(1-1 000)kcV 



ment on board SWIFT.'' In the second group the INTEGRAL/JEM- 
X* and SAXAVFC' rough sensitivities are given compared to the 
HETE 2AVXM* one. Because we concentrate on determining the rel- 
ative number of detections for several missions/instruments (see Sect. 
|4j we will not perform a dedicated calculation of the individual in- 
strumental sensitivities; that goes beyond the scope of our study (see 
Band l2003i for a detailed analysis). 

In order to compare the sensitivities of different missions Band 
<2003 t made emphasis on the need of expressing them in terms 
of a common energy band. Band 1 2003 1 determined the peak pho- 
ton flux threshold (hereafter named Pins, measured in ph cm~^ 
s"^), for several detectors, including HETE 2AVXM, SAX-WFC and 
SWIFT/BAT. INTEGRAL/IBIS and INTEGRAL/JEM-X were not 
considered in that study. In Band l'2003'l all the detection thresholds 
are normalised to the 1 — 1 000 keV energy band, making several as- 
sumptions about the detectors' properties (accumulation time, detector 
response matrix, background model) and GRBs' spectra.' 

Based on Figs. 5 and 7 in Band <2003> we assumed detection 
4 ph cm-2 s"^ for HETE 2/WXM 

s~^ for SWIFT/BAT. These two 
thresholds are valid for a broad range of Ep values of around several 
hundreds of keV, (a typical range for GRBs) and refer to the 1 — 1 000 
keV energy band. We note that Ep is the spectral maximum in the 
vFi, (or equivalent N{E)E^) representation. Many studies have re- 
ported results in the 50 — 300 keV band, mainly because this band 
corresponded to the sensitivity range of BATSE^. Thus, in the present 
work we decided to carry out the sensitivity comparisons in a common 
50 — 300 keV reference energy band. The election of another energy 
band is arbitrary and it would not affect our results, since it would 
only introduce a constant multiplicative factor on the Pins values cal- 
culated for each mission-instrument (as well as on the luminosity, L, 
as it will be shown in equation 2). 

Setting ai ~ a2 ~ 1.5 (see Sect. |3}, we calculated a 

flux ratio p(50-300)kcV/p(l-1000)kcV ^ q qq Poj. ^ pQ^gj. 

law spectrum this detection threshold ratio is independent of 
Ep (which would correspond to a constant straight line in 
Fig. 1 of Band 2003 1. By transforming the assumed 1 — 1 000 
keV thresholds, the derived 50 — 300 keV sensitivities for 



thresholds of Pjj^^^ ^^^^m - 
™d -PswiFT/BAT = 2 ph cm 



HETE 2/WXM and SWIFT/BAT become R 



(50-300)koV _ 



HETE 2/WXM 



Launch scheduled for the second half of 2004. Hereafter, when 
we refer to SWIFT we will focus only on BAT (which provides the 
GRB triggers for the mission) and not on the X-ray and [/V'/Optical 
Telescopes on board SWIFT. See Gehrels et al. <2004> for further de- 
tails. 

■* Joint European Monitor for X-rays, on board INTEGRAL. 
^ Wide Field Cameras on board BeppoSAX. 

* The Wide X-ray Monitor is the most sensitive instrument on 
board HETE 2. The rest of the HETE 2 instrumentation will not be 
considered. 

' In Band <2003t the GRB spectra are described by the Band func- 
tion (Band et aL ll993t , where the photon number flux is given by: 



N{E) = 



I Xnf )"2-'^i( — E — . 

I. ^'"VloOfceV'' VloOfceV. 



— CK2 



exp(Qi — 02' 



E <Ei, 
E> Et 



where Et, — (02 — ai)Eo, and Eo determines the maximum of 
N{E)E^ given by Ep — (2 — qi)£'o. Qi and 02 are the photon 
spectral indexes at frequencies below and above Eb- 

The Burst And Transient Source Experiment; it operated on board 
the Compton Gamma-Ray Observatory between 1991 and 2000. 



0.35 ph cm^^ s^^ and -PgwiFT'/BAT ~ ^-^"^ cm^^ s^^, 
respectively. We note that those sensitivity values refer to the 
peak flux. 

Estimates carried out for the IBAS^ sensitivity yield a 
threshold of - 0.14 - 0.22 ph cm-^ g-i for INTEGRAL/IBIS 
(the ISGRI part) in the 20 - 200 keV energy band 
(Mereghetti et al. "2003'l. Hereafter we assume a conservative 
INTEGRAL/IBIS sensitivity limit of P/^Jj^g^p^^L/iBis = 0.22 
ph cm^^ s^^, which corresponds to ^'inxegral/ibis ~ ^-^"^ 
ph cm^^ s^^, assuming ai = 02 — 1.5. Therefore, the 
INTEGRAL/IBIS capabilities to detect high redshift GRBs rel- 
ative to SWIFT/BAT (presented in Sect. 0} have to be consid- 
ered as a pessimistic estimate of the actual INTEGRAL/IBIS 
potential. 

An alternative, simplified way to verify the photon 
peak flux thresholds for the instruments not considered by 
Band (2003) (like INTEGRAL/IBIS and INTEGRAL/JEM- 
X) would be to assume that the sensitivity threshold Pins for 
a GRB is proportional to the square-root of the background 
count rate and inversely proportional to the square-root of the 
detector area. Thus we have Pins oc y/fl/A, where 57 is 
the sky coverage and A is the detector area. This assumption 
would yield a sensitivity ratio of 3.0 between INTEGRAL/IBIS 
and SWIFT/BAT (see Table [TJ and, therefore, a threshold of 

p(50-300)kcv _ 06 Dh cm-2 s-1 
Mntegral/ibis ~ pn cm s . 

In principle, the most reliable comparison that can be 
done in Table[T]is between SWIFT/BAT and INTEGRAL/IBIS 
(the ISGRI part), because they are based on very sim- 
ilar detector technologies (CdZnTe in SWIFT/BAT, CdTe 
in INTEGRAL/IBIS) and share a similar energy band. 
This simple estimate is a factor of two lower than the 



^iNTECjRAiviBis ~ ^'^'^ cm^^ s^^ threshold assumed for 
INTEGRAL/IBIS. Hence, it supports the use of this threshold 
as a conservative upper limit of the real INTEGRAL/IBIS sen- 
sitivity. 

From the above estimates it is evident that 
INTEGRAL/IBIS will be the most sensitive GRB detec- 
tor (at least in the 15 — 150 keV energy band) ever flown and 
not likely to be matched, sensitivity wise, by any other mission 
within the coming decade. 

Rescaling the HETE 2/WXM threshold with the P,ns oc 
^yn/Amcipe we obtain a sensitivity of -P/ntegral/jem-x ~ 
0.02 ph cm^^s"^. However, the difference in detec- 
tor technologies of INTEGRAL/JEM-X with respect to 
HETE 2/WXM makes this number an uncertain estimate of 
the actual INTEGRAL/JEM-X sensitivity. Furthermore, the on 
flight INTEGRAL/JEM-X performances have been changed 
during the first months of the INTEGRAL mission, so its real 
sensitivity is well above 0.02 ph cm^'^s^^. Besides, the re- 
duced number of bursts that INTEGRAL/JEM-X will detect (a 
very few a year) does not support performing a specific calcu- 
lation aimed at studying its capabilities for high redshift bursts. 

On the other hand an estimate based on Figs. 3 and 5 of 
Band ( 2003 1 yields a relative sensitivity of ^ 3 — 4 between 
HETE 2-WXM and SAX-WFC. This is in agreement with 
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Mission/Instrument 


Area 


Coverage 


Energy band 


Orbit 


Relative 




cm^ 


% of 47r str. 


(keV) 


efficiency 


sensitivity^ 


SWIFT/BAT 


5200 


14* 


15 to 150 


0.6 


1.0 


INTEGRAL/IBIS (ISGRI) 


3000 


0.9* 


15 to 150 


0.8 


3.0 


HETE 2AVXM 


360 


13 


2 to 25 


0.5 


1.0 


INTEGRAL/JEM-X 


1000 


0.1 


2 to 35 


0.8 


19.0 


SAXAVFC 


530 (X 2) 


2(x2) 


2 to 30 


0.5 


3.1 



* Half coded field of view (FOV). 

f Estimate based on the Pms ~ \/Tl/A recipe. 



the HETE 2/WXM vs SAX AVFC relative sensitivity estimate 
given by the Pi„s oc ^/^l/A expression (see TablefO- 

3. Detectability of a high redshift population of 
GRBs 

We have selected the most sensitive missions-instruments 
(present: INTEGRAL/IBIS, HETE 2/WXM, and future: 
SWIFT/BAT) to calculate their capabilities of detecting a high 
redshift population of bursts. To estimate the number of GRBs 
that these missions-instruments will detect we assume that: 

- GRB spectra can be described by power laws (F^ v^"'). 
After exploring the impact of a on the final results (which 
do not change qualitatively) and for simplicity, we have as- 
sumed a value of a = 1.5, typical of GRBs (van Paradijs et 
al.pOOl. 

- The GRB rate is proportional to the star formation rate 
(SFR) in the Universe. The SFR considered is the one given 
by Rowan-Robinson (1999 2001) for z < 5 and the one 



THE STAR FORMATION RATE IN THE UNIVERSE 



calculated by Gnedin & Ostriker ( I1997> for z > 5 (see 
Fig.ID. 

The GRB peak (isotropic) photon luminosity function is 
given by: 



S{L) = 



in 

Otherwise 



being L the peak photon luminosity and /3 the luminosity 
function index. Lmin, Lmax determine the width of the lu- 
minosity function. We have assumed a value of = — 1. 
- Although the effect of several Universe models has been 
checked, we choose to use the most popular cosmological 
parameters in this paper: il„j ~ 0.3, JIa = 0.7, Ho — 65 
km s^^ Mpc^^. 

For the above assumptions the differential GRB detection 
rate at a given peak photon flux P in a detector (ph cm^^s^^) 



is given by the following convolution integral: 

Ngr.b{P) = C^lej^ RGR.BS{L)dL , 



(1) 



where e is the efficiency of the orbit, fl is the instrumental 
coverage of the sky and Rgrb is the GRB detection rate if 
they were standard candles, i.e., Rgrb ^ sfr{z) dv(z) dz^ 
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Fig. 1. The plot shows the SFR in the Universe as a function of 
the redshift. The dashed line represents the SFR derived from 
numerical simulations for z > 5 (Gnedin & Ostriker [1997> . 
The solid line shows the SFR at the z < 5 region based on ob- 
servational estimates (Rowan-Robinson 1999. .2001 J . The tran- 
sition between the two regions has been smoothed. 

being V the comoving volume. The value of the proportion- 
ality constant C is unknown. Fig.|2]displays Ngrb{P) as well 
as the detection thresholds of several high-energy missions- 
instruments. 

The relationship between L, z and P is given by the fol- 
lowing expression: 



P 



L 



47rD(z)2(l + zY 



(2) 



(H 



where D{z) is the comoving distance. In our calculations dif- 
ferent values of a, Lmin, ^max, P, and H.^ are considered. 
Variations of a, (3, Q,^ and Vim do not change the final result 
qualitatively. On the other hand, the variations of Lmin, Lmax 
are more relevant for determining the number of high redshift 
GRB detections. We consider the pessimistic case of a rela- 
tively faint and narrow luminosity function S{L) defined by 
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DIFFERENTIAL NUMHEK OF DETECTIONS (a=l,5, ,!? = -1.0, a„,=0.3, n,,= 0.7) 




RELATIVE NUMBER OF DETECTIONS vs REDSHIFT (ii=i 

1 .0 \ 



S = -1.0, 0^ = 0.3, n,= o.7) 



-3-2-10 1 2 

Log P (ph cm"' s"') 

Fig. 2. Differential peak photon flux distribution of GRBs. The 
soHd curve shows the differential peak photon flux distribu- 
tion if all redshifts are considered, i.e., Ncrb{P)- The dashed 
curves represent the differential peak photon flux distribution 
of GRBs when only GRBs with z > z^dge are taken into ac- 
count, i.e., Ngrb{P, Zedge)- The intention of this figure is to 
show how the Gaussian like solid curve displaces leftwards 
(low P values) as z^dge increases. The dashed lines plotted 
in the figure are shown for the purpose of the example, using 
arbitrary values of Zedge, not marked for simplicity. As Zedge 
increases the closer GRBs (and therefore larger P bursts) are 
left out (right tail of the Gaussian like distribution). The ver- 
tical lines represent the detection thresholds for the different 
missions-instruments, showing the arrows the detectability re- 
gion. 



rmn = 'i^O ph s and L„iax = '^O ph s \ The L,„„ 



L 

and Lmax values used in our study correspond to the narrowest 
S{L) among the ones used by Lamb & Reichart (2000 1, and 
consistent with the observed peak photon luminosity distribu- 
tion (the measured S{L) is at least 1.7 times wider; Stern et 
al. I2002> ^°. The assumption of a wider S{L) would broaden 
the Gaussian like curves (both solid and dashed) displayed in 
Fig. 121 extending their tails to higher P values and therefore 
enhancing the number of GRBs detected at very high redshift. 

We can calculate the contribution to (1) by GRBs with red- 
shift larger than z^dge (see dashed curves of Fig.|2j, using: 



GRB 



poo 

{P, Zedge) =Cne H{z{L) )RGB.BS{L)dL , 

Jo 

where H{z{L), z^dge) is a step function that vanishes unless 
z{L) > z,dge- Obviou sly, Ngrb{P) = ^GflB(-P,0), and 

Our Lmin and Lmax values are also in agreement with the em- 
pirical peak photon luminosity function used as reference by Lamb & 
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Fig. 3. Relative number of detections as a function of the 
redshift. This plot shows, for several high-energy missions- 
instruments, the fraction of the detected GRBs that have a red- 
shift larger than Zedge- 



< 1. Last, we can calculate the number of 



Ngrb{P) 

GRBs detected above a given instrumental peak photon flux 
threshold Pins that have redshifts larger than Zedge- 



GRB 



{Pins, Zedge) ^ J NgRb{P, Zedge)dP . 



Reichart f200OT; log L = 58.1 ± 0.7 



Not knowing the proportionality constant C of (1), we can 
not derive an absolute value for NGRB{Pins, Zedge)- However, 
we can determine the relative quantity ^'no'j^bI'p- '^"q)'^ , which 
provides us the proportion of detections that have a redshift 
larger than Zedge (see Fig.|3. 



4. INTEGRAL/IBIS vs SWIFT/BAT; comparison of 
the number of GRB detections 

Fig. |3] shows the relative number of detections (given 
by ^GHg(^""-^-'^g°) ) as a function of the redshift, for 
HETE 2/WXM, SWIFT/BAT and INTEGRAL/IBIS. As it is 
shown, - 15% of the GRBs detected by INTEGRAL/IBIS 
will have a redshift larger than 4. For SWIFT/BAT the z > 4 
population will be ~ 10% of the total number of detections. 
In the case of INTEGRAL/IBIS the tail of NoRB{P.,^s,z.dg.) 

NGRB{Pins,0) 

extends even up to redshifts of Zedge > 7. SWIFT/BAT and 
HETE 2AVXM will detect a closer population of bursts, spe- 
cially HETE 2/WXM. HETE 2/WXM is the less sensitive 
GRB mission, being constrained to detect bursts with red- 
shifts z < 4.3. This prediction is in agreement with the maxi- 
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mum redshift measured for a HETE 2/WXM GRB. ' ' Therefore 
we will not consider HETE 2/WXM in the following study, 
aimed at calculating the relative number of GRB detections 
as a function of the redshift. We will concentrate on compar- 
ing INTEGRAL/IBIS and SWIFT/BAT capabiUties. Besides, 
as we previously noted, the similar energy range and detector 
technologies of INTEGRAL/IBIS and SWIFT/BAT suggest a 
reliable comparison. 

For determining the relative number of detections between 
two missions, A and B, the following calculation has to be per- 
formed: 

, . _ ^a^aJp^ Jo°° H(z(L),z^i,^)RGB.BS{L)dLdP 

jA-B[Zedge) - nseB H{z{L),z,a,,)RGRBS{L]dLdP ' ^^^^ 

This function will give the relative number of GRB 
detections with z > z^dge- We have applied (2) to 
derive the detection ratio between INTEGRAL/IBIS and 
SWIFT/BAT as a function of the GRB population redshift 

(/lNTEGRAL/IBIS-SWIFT/BAT(^edge)) ■ 

As it is shown in Fig. |3 for Zf,dge < 5.6 
/iNTEGRAL/iBis-swiFT/BAT < 1- In Other words, at low 
redshifts the large field of view (FOV) of SWIFT/BAT, 
in comparison to INTEGRAL/IBIS, dominates the num- 
ber of detections. On the other hand, for Zedge > 5.6, 
INTEGRAL/IBIS sensitivity becomes the dominating factor 
and /iNTEGRAL/iBis-swiFT/BAT > 1- Thus, from the point 
of view of the ground based strategy of detecting z > 6 after- 
glows, it might be more efficient to observe INTEGRAL/IBIS 
GRBs than to observe SWIFT/BAT bursts. If we are interested 
in using GRBs to study the reionisation epoch that occurred at 
2 ~ 7 (Loeb & Barkana r2001> . then it would be relevant to pri- 
oritise the foUow up of INTEGRAL/IBIS GRBs. Nonetheless 
we emphasise that INTEGRAL/IBIS high redshift detectabil- 
ity predictions are based on low number statistics, so they are 
subjected to large fluctuations. 

The self consistency of our procedure can be checked 
by studying the prediction of expression (3) for the partic- 
ular case when Zedge = 0. If we consider z^dge = 0, 
then /integral/ibis-swift/bat(0) gives us the fraction of 
GRBs detected with z > 0, i.e., considering all the detec- 
tions independently of their redshifts. We predict a value of 
/integral/ibis-swift/bat(0) = 1/9.7 for the fraction of 
the total number of GRBs detected by INTEGRAL/IBIS com- 
pared to SWIFT/BAT (see Fig. Hat z^dge = 0). The number of 
GRBs cun-ently detected by INTEGRAL/IBIS is - 1 GRB per 
month (Mereghetti et al. i2003j . whereas the last updated esti- 
mates for the SWIFT/BAT detection rate yield - 110 GRBs 
per year (Hevl l2003t . This gives a fraction of ^ 1/9.2 for the 
number of detected GRBs, in agreement with our prediction. 

Although INTEGRAL/JEM-X FOV and sensitivity are 
less suitable than the ones of INTEGRAL/IBIS to detect 
GRBs, the spectral peak of the high redshift GRBs (usu- 
ally at 100 — 1200 keV) will be in the detection range of 

" Among the 12 HETE 2AVXM GRBs with confirmed spec- 
troscopic redshifts to date, the maximum redshift is reached for 
GRB 030323 at 2 = 3.372 (Vreeswijk et al. l2004l . 



RELATIVE NUMBER OF DETECTIONS BETWEEN 
INTEGRAL/IBIS AND SWIFT/BAT 
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Fig. 4. The plot shows the GRB detection ratio of 
INTEGRAL/IBIS with respect to SWIFT/BAT as a function 
of the GRB population redshift. /integral/ibis-swift/bat 
increases with redshift because of the better sensitivity for 
high redshift GRBs of INTEGRAL/IBIS. At redshifts z < 
5.6 the larger field of view of SWIFT dominates, making 
/integral/ibis-swift/bat < 1- 

INTEGRAL/JEM-X. So given that INTEGRAL/JEM-X is co- 
aligned with INTEGRAL/IBIS, it might be also useful to detect 
the redshifted prompt 7-ray component. 

5. NIR observations for INTEGRAL GRBs 

The most significant inconvenience to follow up INTEGRAL 
GRB afterglows is the large fraction of the mission's time de- 
voted to scan the extincted Galactic centre/plane. An additional 
problem (not specific to INTEGRAL like the previous one) 
for detecting high redshift GRBs comes from the Ly-a blan- 
keting effect that strongly attenuates the radiation observed at 
A < 1216 X (1 + z) A. Both drawbacks can be mitigated if 
ground based searches are carried out in the NIR, specially in 
the ii'-band, where i) the Galactic extinction drops by a fac- 
tor ^ 6 compared to the optical, and ii) the Ly-a blanketing 
problem is eliminated up to very high redshifts (z ^ 17). 

There are further advantages to searching for high redshift 
afterglows in the NIR. Because of their fading behaviour (typi- 
cally, GRB afterglows decay following a power law of the form 
Fi, ^ t^*, where t is the time since the onset of the 7-ray event 
and S the decay index) the time dilation effect increases the ob- 
served flux at a fixed time of observation after the GRB. Thus 
current and future ground based NIR facilities could be able to 
detect high redshift afterglows that are invisible in the optical 
bands. 
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As an example, in Fig. |5l we display the very well sam- 
pled specti-al energy distribution (SED) of the GRB 000301c 
afterglow (Jensen et al. 1200 U for different redshifts, once the 
Ly-a blanketing absorption has been modeled (Madau'1995 1. 
The photometric point flux densities (F^; see dashed line in 
Fig. |5} have been fitted following an expression of the form 
Fi, ^ v^' X 10^° '*^'^, where rj is the afterglow power law spec- 
tral index and A^, is the absorption in the GRB host galaxy at a 
frequency of ly. has been parameterised in terms of Ay fol- 
lowing a typical SMC extinction law (Pei' 1992V For illustra- 
tion purposes we have overplotted, for different exposure times, 
the JHK-hwA 5(t sensitivities foreseen for the Espectrografo 
Multiobjeto InfraiTojo (EMIR) (Balcells'1998't, currently being 
build for the 10 m Gran Telescopio Canarias (GTC) (Rodriguez 
Espinosaetal. [T998b '^ 

We would like to emphasise several points. First, this sim- 
ple redshifting exercise does not require any assumption on 
the star formation rate, luminosity function, etc. . . Second, the 
GRB 000301c afterglow was not specially bright, so an intense 
afterglow (like GRB 990123; Castro-Tirado et al. 1999) could 
be detectable even at higher redshifts. Third, the photometric 
points displayed in the figure correspond to a fairly late epoch 
observation (they were acquired ~ 1.06 days after the burst, 
time measured in the observer's frame), so afterglows observed 
a few hours (even minutes) after the burst would be reachable 
by smaller facilities (i.e., NIR robotic telescopes like REM; 
Zerbi et al. l200Tl 12003). 

As it is shown in Fig. [Sj rapid NIR observations of after- 
glows performed with 10 m class telescopes (e.g., the GTC 
equipped with EMIR) would be dominated by the Ly-a blan- 
keting effect since, above a given exposure time threshold 
{Texp ~ 900 s for the particular case of GRB 000301c and 
the GTC), the detectability is basically independent of the ex- 
posure time employed in the observations. 

At high redshifts the observer would be samphng the UV 
part of the afterglow synchrotron spectrum where the extinc- 
tion of the GRB environment is expected to be severe. Thus 
one of the mayor problems to predict the NIR detectability of a 
high redshift afterglow would come from the, so far unknown, 
extinction law and dust content describing both the local GRB 
environment and the large scale line of sight properties within 
the host galaxy. Furthermore, the UV opacity in the environ- 
ment close to the GRB progenitor might be time dependent. 
Hence, a detailed quantitative prediction of the number of af- 
terglows detected in the NIR (beyond the scope of this work) 
would require modeling the physical processes describing the 
UV absorption of the afterglow radiation and its evolution at 
different scales (from several AU to kpc) around the GRB pro- 
genitor. 

6. Discussion and conclusions 

In the present paper we study the capabilities of several 
missions-instruments to detect high redshift GRBs. At low 

For additional information, please, visit the following sites: 
http : / /www ■ ucm .es/info/ emir/index_e ■ html| and also: 
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• Photometric points 

♦ Texp = 900s 
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Fig. 5. The plot shows the evolution of the observed flux of 
a typical GRB afterglow (GRB 000301c) when it is redshifted 
from z — 2.0404 (spectroscopic redshift; Jensen et al. 1200 U 
up to z = 22. The circles show the photometric measurements 
of the afterglow from the U to the iiT-band, taken 1.06 days af- 
ter the GRB. The dashed line represents the SED fitted to the 
photometric points (see main text). The rhomboids, triangles 
and squares show the JHK-hwA sensitivities {5a) foreseen 
for the future GTC 10m telescope {+ EMIR instrument), for 
several exposure times. As it can be seen afterglows can be 
detected in the three JHK bands up to z = 9, with only a 
exposure time of 900 s per filter A'-band counterparts would 
be detectable up to z ^ 17. For T^xp > 900 s the detectabil- 
ity of the GRB 000301c afterglow with the GTC (h-EMIR) is 
dominated by the Ly-a blanketing effect, and is basically in- 
dependent of the exposure time employed in the observations. 
The Ly-a blanketing absorption has been modeled following 
Madau ( fT995t . 



and moderately high redshifts (0 < z < 5.6) the large 
FOV of SWIFT/BAT, compared to the instrumentation on 
board INTEGRAL, dominates the number of GRB detections. 
Nonetheless, the better sensitivity of IBIS makes INTEGRAL 
more efficient detecting GRBs beyond z ^ 5.6. 

In the most popular cosmological models, the first sources 
of light began at a redshift of z 30 and reionised most of the 
Universe by z ~ 7 (Loeb & Barkana 2001 1. Thus the study of 
the Ly-a forest, present in the high redshift INTEGRAL/IBIS 
afterglow spectra, might constrain the epoch of such reioni- 
sation (Gunn & Peterson 1965 1, probing the ionisation state 
of the intergalactic medium (IGM) as a function of redshift. 
According to our estimates this kind of studies will be less pro- 
ductive if SWIFT/BAT (even more HETE 2/WXM), rather than 
INTEGRAL/IBIS, GRBs are followed up. 

Further, optical/NIR spectroscopy of high redshift after- 
glows would provide additional information on the primitive 
IGM and ISM. High resolution spectroscopy of high redshift 
afterglows might reveal metal absorption lines, which could 
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trace the metal enrichment history of the Universe. In princi- 
ple this is the same technique as the one used for studying the 
damped Ly-a systems (DLAs), which are located in the line 
of sight of quasars. However GRBs are cleaner probes than 
quasars, because a brief phenomenon such as a GRB would 
not modify the state of the environment at large distances as 
quasars do (the Stromgren sphere of high redshift quasars is 
of the order of Mpc; White et al. 2003j. GRBs would allow 
to study the ISM of basically unaltered host galaxies, highly 
ionised if they had harboured a quasar during at least one mil- 
lion years. Moreover, GRBs are brighter (although for a short 
time) than quasars, so GRBs can probe the ISM of dustier 
galaxies. 

Theoretical studies show that the birthrate of Population III 
stars produces a peak in the SFR in the Universe at redshifts 
16 < -z < 20, while the birthrate of Population II stars pro- 
duces a much larger and broader peak at redshifts 2 < z < 10 
(Valageas & Silk l 19991 . If GRBs are produced by the collapse 
of massive stars, they are expected to occur at least at z 10 
and possibly up to z ~ 15 — 20. The mere detection of very 
high redshift GRBs would give us for the first time direct con- 
firmation of the existence of the earliest stellar generations. 

Rapid NIR observations are a complementary strategy to 
study high redshift INTEGRAL/IBIS GRBs. iC-band observa- 
tions performed with current and planned ground based NIR fa- 
cilities might detect GRBs up to the theoretical limit imposed 
by the Ly-a blanketing (z ^ 17). Realistically however the 
z ~ 17 upper limit ought to be decreased due to possible UV 
absorption present in the GRB host galaxy. 

In conclusion, the INTEGRAL/IBIS 's capabihties for 
studying GRBs are complementary to the ones of missions- 
instruments like SWIFT/BAT and HETE 2/WXM, spe- 
cially devoted to prompt localizations of GRBs. Whereas 
SWIFT/BAT and HETE 2/WXM would detect more GRBs 
than INTEGRAL/IBIS, the latter might detect high redshift 
GRBs unreachable to the earlier. Future works might invert 
the logical outflow followed in the present paper Once spec- 
troscopic redshifts have been measured for a large sample of 
GRBs (and therefore the number of detected GRBs as a func- 
tion of redshift is known), the equations could be inverted in 
order to obtain the SFR evolution with redshift. INTEGRAL, 
and specially IBIS, could be a very valuable tool to trace the 
SFR rate in the early Universe. 
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